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In a recent paper |l| we considered the radiative decay 
4> — » 70,0 / /o in the molecular (KK) model of the scalar 
mesons (ao(980), /o(980)). In particular, we showed that 
there was no considerable suppression of the decay am- 
plitude due to the molecular nature of the scalar mesons. 
In addition, as a more general result we demonstrated 
that, as soon as the vertex function of the scalar me- 
son is treated properly, the corresponding loop integrals 
become very similar to those for point-like (quarkonia) 
scalar mesons, provided reasonable values are chosen for 
the range of the interaction. We also confirmed the range 
of order of 10 -3 -j- 10 -4 for the branching ratio obtained 
in Refs. @,i,H within the molecular model. 

As a reaction to our work a paper appeared [B| , where 
the authors criticize our results and claim that our paper 
[H is "misleading" . Specifically, they dispute our findings 
that the transition amplitude 4> — > K + K~ — > 7a //o is 
governed by low kaon-momenta (nonrelativistic kaons) 
in the loop. In order to support this conjecture they 
present numerical results that supposedly demonstrate 
that "ultrarelativistic kaons determine the real part of 
the <p ~ * K + K~ — ► 7ao//o amplitude". The dominance 
of such contributions of "kaon high virtualities" is then 
interpreted as support for a compact four-quark nature 
of the scalar mesons. 

In this comment we want to point out a fundamen- 
tal flaw in the calculations presented in Ref. [B| which, 
in turn, invalidates the criticism raised in that pa- 
per. Namely, in order to demonstrate that the high- 
momentum components determine the <f> — > K + K~ — > 
7 a o / fo amplitude the authors of [H introduce a momen- 
tum cut-off in the relevant integrals. However, in doing 
so gauge-invariance gets violated. As will be shown be- 
low, large momentum contributions appear only in this 
induced gauge-invariance-violating term and are there- 
fore of no physical significance. 

To keep our argument self contained we briefly repeat 
the essentials of the formalism. As a consequence of 
gauge-invariance the full matrix element for the <f> — » 7 5* 
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FIG. 1: Diagrams contributing to the amplitude of the radia- 
tive decay </> — > yao/fo- 

(S = ciq or fo) decay, M v , can be written as 

M v = eg0g5 /(my, m s )[e v (p ■ q) - p u (q ■ e)] 

= eg^gse^J^v, (1) 

where p and q are the momenta of the <p meson and the 
photon, respectively, m is the kaon mass, g<f, and gs are 
the 4>K + K~ and SK + K~ coupling constants, and e v is 
the polarization four-vector of the 4> meson. The masses 
of the 4> meson and the scalar are denoted by my and 
ms, respectively. The function J(my,ms) has a smooth 
limit for q — > 0. As a consequence of gauge-invariance 
the amplitude |T]) is transverse, M v q v — 0, and is pro- 
portional to the photon momentum; especially it vanishes 
for q — > 0. The form {l} is well known. Details can be 
found, for example, in Refs. [1, 0, BHEl- 

For point-like scalars, only diagrams (a)-(c) of Fig. Q] 
contribute. If the scalars are regarded as extended ob- 
jects, a vertex function needs to be introduced at the 



2 



KKS vertex. Then gauge invariance demands the inclu- 
sion of a diagram of type (d). For general kinematics a 
proper construction of this additional term is quite in- 
volved (see Ref. [H| where we list a few of the papers 
devoted to this subject) and contains some ambiguity. 
However, for soft photons, all the different recipes give 
the same result up to corrections of order (q/f3) 2 that will 
be dropped. Here 1 /f3 denotes the range of forces — for 
the case of interest one may use [3 ~ m p , where m p de- 
notes the mass of the lightest exchange particle allowed, 
namely that of the p meson We may then as well 
use the method suggested in Ref. [l(| that is based on 
minimal substitution considerations. For more details on 
the issue of gauge invariance for the reaction considered 
here, see Ref. [12j . 

After this introduction let us discuss the main formula 
of our paper It is argued there (and confirmed by 
actual calculations) that, since the amplitude is finite 
even for the point-like limit, the range of convergence of 
the integrals involved is defined only by the kinematics 
of the problem. In particular, if both masses, i.e. that 
of the vector and of the scalar meson, are close to the 
KK threshold, the integrals converge at fco ~ m, and 
thus for non-relativistic values of the three-dimensional 
loop momentum k, \k\ <C m. This allows us to perform 
a nonrelativistic reduction of the amplitude in the rest- 
frame of the (f> meson. The integrals in question for the 
individual graphs of Fig. 1 are (note that jj^ = J^)'- 
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where Ey = my — 2m, Es — ms — 2m, and K = k — 
^q. The last integral can be rewritten by performing an 
integration by parts: 
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Here, contrary to Ref. with the last term we kept 
the surface integral that emerges in the calculation. In 
order to investigate the range of momenta relevant for 
the loop integrals in Ref. H, a momentum cut-off A was 
introduced. We follow this prescription and write the full 
transition current as: 



Jifc(A) = Jifc(A) +S ik R(A) 



where 
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contains the above mentioned surface term and 
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For later convenience we used energy conservation to 
replace Es via Es = Ey—q 1 . For A — >• oo Ji k (A) matches 
to the formula used in Ref. [l| to calculate the matrix el- 
ement for (j> — > 7ao//o- We checked that this sum of inte- 
grals converges for non-relativistic kaon momenta. This 



1 Contrary to the claim made in reference [5] of Ref. [e| energy 
and momentum conservation are maintained in the calculations of 
Ref. Q. 



finding was confirmed in Ref. Q. 

For illustration we choose a particular form of r(k), 
namely T(fc) = (3 2 /{k 2 + (3 2 ), and study that part of J ik 
proportional to the structure Si k (according to Eq. |T]) ex- 
actly this structure contributes to the decay amplitude in 
the 0-meson rest frame). In Fig. [2] we plot the behaviour 
of the integrand j(k) ( Im{Ji k ) — 5i k L j(k) dk), as 
a function of k (note, that the integrand in the simi- 
lar integral for Re(Jifc) contains S(k 2 — mEy), such that 
k = y/mEy w 0.12 GeV). From Fig. [5] one can see that 
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FIG. 2: The behaviour of the integrand j(k) of Im( Jik), as 
a function of the kaon momentum floating in the loop. We 
chose m s = 0.98 GeV, m v = 1.02 GeV, m = 0.495 GeV, and 
four values of the parameter f3: (5 — 0.4 GeV (dotted line), 
(3 = 0.6 GeV (dashed line), = 0.8 GeV (thin solid line), and 
f3 = oo (thick solid line). 



the integral indeed converges at non-relativistic values 
of the kaon momentum, regardless of the value of the 
finite-range parameter /3 — the latter plays no role for 
the convergence. 

In Ref. [l[ the last term of Eq. (g]), R(A), was dropped, 
for it vanishes exactly for A — > oo 2 . In Ref. Q, however, 
it is argued that this term should be kept and that it con- 
verges only for very large values of A, which means that 
the corresponding integral acquired contributions from 
very large momenta. The contribution of those large mo- 
mentum components is then taken as a proof that only 
if the scalars are very compact objects, a sizable contri- 
bution from the loop can emerge. Notice that, even for 
finite values of A, J,fc(A) vanishes for q — > 0, as required 
by the general structure given in Eq. ([1]) . However, since 
i?(A) is independent of the photon momentum q, it gives 
a nonvanishing contribution to Jik even for q = for all 
finite values of A. Therefore this term violates gauge- 
invariance. Thus, by introducing a sharp cut-off into the 
problem the authors of Ref. [5J produced a term that 
violates gauge-invariance 3 . Since the whole argument 
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true we only need to demand that 



3 With sharp cut-off, gauge-invariance of the amplitude can be 
restored by a subtraction at q = 0. Obviously, this procedure is 
equivalent to omission of the last, (/-independent term in Eq. Ht . 



presented in Ref. [B| is based on this term, it bears no 
physical significance. 

We therefore conclude that all results of Ref. [l[ are 
valid. In particular, there is no strong suppression of kaon 
loops by the scalar wave function. Regardless of this, it 
should be stressed that the data for (f> — > 7a //o [13|] is 
very sensitive to the nature of the light scalar mesons, for 
it allows direct access to the effective coupling constant 
g c g of the scalar to the kaons. As was shown in Ref. Q 
this coupling is a direct measure of the molecular contri- 
bution of the scalar mesons. 

As a reply to this comment, the authors of the com- 
mented paper [B[ state that, with a properly regularised 
amplitude T„ M (up to the overall normalisation T Vfl coin- 
sides with our J„ M introduced in Eq. {T]) above), "the 
regulator field contribution is caused fully by high mo- 
menta {M — > oo) and teaches us how to allow for high 
K virtualities in gauge invariant way" 1 1511 . To exem- 
plify their statement, the authors of Ref. [15j employ the 
Pauli-Villars regularisation in the full relativistic expres- 
sion for the amplitude, 

T{cp(p) -> 7 {a (q)/fo(q)],M} = {^{^T^ip, q), 

(7) 

where the overline marks the regularised quantities. In 
particular, for a quantity 
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0= / d*rf(r,p;m), 



the corresponding regularised integral reads: 

0(M) = f d*r[f(r,p;m) - f(r,p;M)}, (9) 



with A/ being the regulator mass; the physical amplitude 
corresponds to the limit M — > oo. Then, by an explicit 
calculation, the authors of Ref. 15] find that 



e"(0K(7)lV (P,P) 

= e<W( 7 ) {T™(p,p)-T%(p,p)) 

= (e(0)e( 7 ))(l - 1) = 0, (10) 

and thus they conclude that the regularised amplitude is 
gauge invariant, and this is due to the final subraction 
coming from the regulator piece. To have a deeper in- 
sight, consider T^(p,p) and extract the coefficient of the 
structure g^ v . Then, after using the Feynman parameti- 
sation. one has: 
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where C is an unimportant numerical coefficient and 
R 2 = z(l — z)p 2 — m 2 . The integrand can be rewritten 
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in the form: 
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It is tempting now to perform averaging over the angu- 
lar variables first, substituting r^rv — > -gT g^v Then, 
naively, the first term in the integrand (fl~2|) vanishes, 
whereas the remaining, second term gives a finite con- 
stant independent of the mass m. Conclusions made by 
the authors of Ref. 15j are based on this result, and, fi- 



nally, they notice that "the finiteness of the subtraction 
constant hides its high momentum origin" . However, the 
analysis performed above has a flaw. Indeed, although 
the angular integration of the first term in (fl"2"|) gives 
zero, the remaining radial integral diverges. Therefore, 
one deals with an undefined expression of the kind x oo. 
In order to resolve this issue, it is important to deal with 
finite integrals. To this end, we evaluate the integral 
in d = 4 — e dimensions. The radial integral is finite 
now, whereas the angular integral gives the substitution 
TyXv — > -TZTJ^dtiv- It is easy to check that, after taking 
the limit e — > 0, the contribution of the first term does 
not vanish any more but, on the contrary, it cancels the 
contribution of the second term in (|12p. Thus 
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and T%(p,p) = 
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individually. Thus, the apparent contribution of high 
kaon virtualities found in Ref. [151 ] appears as a result 
of misusing of the Pauli-Villars regularisation scheme, 
when a regularised finite integral is artificially split into 
divergent parts and each part is considered separately. 
In other words, the physical amplitude is given by a sum 
of several divergent integrals and, in any correct regu- 
larisation scheme, these infinite parts cancel each other 
exactly. 

Finally, had the infinitely high momenta been relevant, 
as advocated in Refs. [5j, |l5(, then the pointlike limit 
would have rendered useless, as there is no infinitely com- 
pact states in nature. This becomes clear when a scalar 
form factor is introduced, in which case neither dimen- 
sional nor Pauli-Villars regularisation is needed. In order 
to maintain gauge invariance one is forced either to in- 
clude the diagram (d) of Fig. 1, or to subtract explicitly 
the amplitude, given by the sum of the diagrams (a)- 
(c), at q — ► 0, with both procedures being equivalent. 
Notice that gauge invariance requires the subtraction of 
the original amplitude, and not the one with mass m re- 
placed by the (infinitely large) regulator mass M. The 
latter observation invalidates the claim made in Ref. 15]. 
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